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We present a Gaussian ensemble of random cyclic matrices on the real field and study their spectral 
fluctuations. These cyclic matrices are shown to be pseudo-symmetric with respect to generalized 
parity. We calculate the joint probability distribution function of eigenvalues and the spacing dis- 
tributions analytically and numerically. For small spacings, the level spacing distribution exhibits 
either a Gaussian or a linear form. Furthermore, for the general case of two arbitrary complex 
eigenvalues, leaving out the spacings among real eigenvalues, and, among complex conjugate pairs, 
we find that the spacing distribution agrees completely with the Wigner distribution for Poisson 
process on a plane. The cyclic matrices occur in a wide variety of physical situations, including 
disordered linear atomic chains and Ising model in two dimensions. These exact results are also 
relevant to two-dimensional statistical mechanics and //-parametrized quantum chromodynamics. 
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With thepseudo-Hermitian extension of quantum me- 
chanics [l|-|3 1 it has become possible to develop a number 
of new ideas, opening thereby interesting and important 
directions of investigation. One of these advances has 
been in random matrix theory where pseudo-unitarily 
invariant ensembles were presented Q that exhibit com- 
pletely different kind of level repulsion as compared to the 
ensembles known Thus, physical systems that vio- 

late parity and time-reversal invariance (T^T-symmetric) 
exhibit level repulsion that could be linear or ^ — 5 log S 
where S is the nearest-neighbour spacing of levels. How- 
ever, an explicit analysis has been done only for an en- 
semble of 2 X 2 matrices. 

Irr this Letter, we present random matrix theory 
(RMT) oi N X N cyclic matrices with real elements. 
As we shall show, these matrices are pseudo-symmetric 
with respect to "generalized parity" . Such matrices arise 
in very significant contexts, the celebrated example be- 
ing that of Onsagcr solution of two-dimensional Ising 
model They arc encountered in the treatment 

of linear atomic chains with Born-von Karman bound- 
ary condition [l^l and in understanding overlap matrices 
for molecules like benzene. These matrices also occur 
as transfer matrices in the theory of disordered chains 
[m and in the general context of wave propagation in 
one-dimensional structures [T^]. In the latter example, 
generally, matrices of second order occur - thus, our ear- 
lier results [13[ throw light on the fluctuation properties 
of the eigenvalues. Cyclic matrices also appear in the 
context of phase transitions in the spherical model p^ . 
In all these varied instances, as soon as there is a random 
parameter (e.g. external field or a random coupling in the 
example of Ising model) , the level correlations dictate the 
long time tails of the time correlation functions which, in 



turn, relate to the relaxation of these systems when they 
are perturbed from thermodynamic equilibrium [l^. 

RMT appears in seemingly unrelated problems in 
physics and mathematics ranging from growth models, 
directed poly mers, random sequences, to Riemann hy- 
pothesis [la - llq . Also, the study of random matrices 
has been related to quantum chaos and exactly solvable 
models in a remarkable way [T8l-l20|. Generically, the 
statistics of spectral fluctuations of classically integrable, 
pseudointegrable, and chaotic systems follow respectively 
the general features of Poisson, short-range Dyson model 
[21I I22I or Semi- Poisson [2^ , and Wigner- Dyson ensem- 
bles [2J|. However, for the physical situations occur- 
ring in two-dimensional statistical mechanics where time- 
reversal and parity are violated [25|, [27l - [29j , there is no 
general understanding of the statistical nature of spectral 
fluctuations [13, HH . Perhaps the first example of a bil- 
liard system with a T'T-symmetric (violating V and T) 
Hamiltonian was a particle enclosed in a rectangular cav- 
ity in the presence of an Aharonov-Bohm flux line (3^ . 
For this classically pseudointegrable system, the spectral 
statistics of quantum energy levels was found to exhibit 
level repulsion that is distictly different from the stan- 
dard RMT 0. For these class of systems, an important 
step was taken in and the present work takes us to 
show the nature of these fluctuations in x iV cyclic 
matrices. The general case oi N x N random pseudo- 
Hermitian matrices remains open, however. 

Let us consider an N x N cyclic matrix with real ele- 
ments, {a;}: 
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It is important to note that this matrix is, in fact, pseudo- 
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Hcrmitian (pseudo-orthogonal) with respect to rj 
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that is, 
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Since t)^ = identity, \, r] is introduced here as "general- 
ized parity" . Thus, we have an ensemble of random cyclic 
matrices (RCM) that is pseudo-orthogonally invariant in 
the sense of ([3]). There arc two distinct scenario with 
respect to time- reversal, T and parity, V: (a) standard 
case where T and V are preserved, this case is trivially 
PT-symmetric, and, (b) the case of PT-symmetry where 
T and V both are broken. In case (a), one may study 
the fluctuations properties of energy levels after classify- 
ing the eigenfunctions according to definite parity (odd 
or even) ; however the case (b) belongs to a different class 
altogether. Whereas case (a) corresponds to the invariant 
ensembles of random matrix theory jlj , case (b) has not 
been fully studied, only some partial results exist d, [l^ 
and RCM belong to this case. To our knowledge, the 
discrete symmetries for operators represented by cyclic 
matrices are clearly spelt out here for the first time. Due 
to this generality, our final results are expected to be rel- 
evant for a wide variety of physical situations occurring 
in anyon physics [25[, i/-parametrized quantum chromo- 
dynamics, fractional quantum hall systems (26j . etc. 
The eigenvalues of M are given by [ssj 



N 



El =^apcxp— (p-l)(?-l); 



(4) 



(l = 1,2,...,N), the maximum real eigenvalue being 
Oi. The diagonalising matrix is given by [l^ 



Lcxp— (j-l)(Z-l). 



(5) 



We consider a Gaussian ensemble of cyclic matrices with 
a distribution, 



P(M) - exp-A tr (M'^M) 



(6) 



where A sets the scale (of energy, for instance). 

For the sake of simplicity, we present the analysis for 
an ensemble of 3 x 3 matrices. We would like to ob- 
tain the joint probability distribution function (JPDF) 
of eigenvalues because all the correlations are related to 
it. Also, we would like to show results on the spacing 
distribution as they enjoy a central place in discussions 
in quantum chaos, universality arguments, and rule the 
dominant long-time tail in correlation functions. We im- 
mediately see that tr M'''M = 3(af + a\ + al). In ef- 
fect, we have P({a,}) = (M) ^ g-^^^- . There are 



three eigenvalues - one real, Ei = J^i '^i ^^d a com- 
plex conjugate pair, (i?2,£'2)- We may define spac- 
ing as 5*23 \E2 — E^l ~ \/3(a3 — 02) as well as 

5i2 \Ei - E2\ = ||(a2 + as) + ^(02 - as)!- Ob- 
viously, S12 = 5i3. The JPDF of eigenvalues P{{Ei\) 
can be written as 



P{E,,E2,E*) 



A 



^~A{Ei+2\E2\^ 



(7) 



With this JPDF, spacing distributions can be found 
[3^ . Spacing distribution for the complex conjugate pair, 
-Pec (5*23) is given by 

3 

PcciS23) = l[da,P{{a,})S{S23 - V3\ as - a2\) 



(8) 



Using this, we may define an average spacing, 6*23 
through the first moment and obtain finally a normalized 
spacing distribution in terms of the variable z = 823,/ 823'- 



Vcc{z) 



(9) 



Similarly, the spacing distribution, Prc{Si2) is ob- 
tained: 



(10) 



Mean spacing turns out to be 512 = fv^*^ where 
c = 2F1 [f , |, 1, i] = 1.31112... Defining z = ^12/^, 
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Prc{z) = c zexp I -y^c z 
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(11) 



We can now make following observations ; (i) the 
Gaussianity of Pcc{z) implies that there is no level re- 
pulsion among the complex conjugate pairs, at the same 
time there is no attraction, there is no tendency of clus- 
tering as in Poissonian spacing distribution; (ii) real and 
complex eigenvalues display linear level repulsion. These 
results are also borne out by the numerical simulations 
in Fig. [Hand Fig. H 

For the general case oi N y. N matrices, we need to in- 
vert . This inversion leads us to the following relation: 



^iiEi 



(12) 



where S^; = ^(*-i)(Af-(i-i)) and uj = e^""'^ is a root 
of unity. S is a symmetric matrix and = Nr]. Em- 
ploying these relations, we can find af , and hence the 
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FIG. 1: Probability distribution of the absolute spacing be- 
tween the complex conjugate pair of eigenvalues of a Gaussian 
ensemble of 3 x 3 cyclic matrices. The numerical result ob- 
tained by considering 10000 realizations agrees with the an- 
alytic result ((9|. The Gaussian spacing distribution may be 
interpreted to give an accumulation of eigenvalues resulting 
in a maximum at zero spacing, but no tendency to cluster as 
the first derivative is zero. This is different from a Poisson 
distribution. 
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FIG. 2: Probability distribution of the absolute spacing be- 
tween a real and a complex eigenvalue of a Gaussian ensemble 
of 3 X 3 cyclic matrices. The numerical result obtained by 
considering 10000 realizations for 3x3 matrices and 1000 re- 
alizations of 100 x 100 matrices agrees with the analytic result 
(|ll|l . We observe a linear level repulsion near zero spacing, 
however the result is distinctly different from the Wigner sur- 
mise for GOE. 



where Ei and En__^_-^ real and the rest of the eigenvalues 
may be complex. For odd N, the above result will hold 
except that there will be only one real eigenvalue, Ei 
and the summation in the second term will extend over 
all i except 1. Employing this general result on JPDF, 
we can now calculate the spacing distributions for the 
general case. There are three cases : (i) spacing among 
the complex conjugate pair of eigenvalues is found to be 
distributed again as a Gaussian; (ii) spacing between a 
real and a complex eigenvalue is distributed according to 
(jlip : (iii) two complex eigenvalues, Ej = xj + iyj and 
Ek = Xk + iyk are spaced according to 

/n, d^Ed'^E,pmm\E, - Ek\ - s) 

Jllid^E,d^E,P{{E,}) ' ^'^^ 

which reduces to the following integral on change of vari- 
ables, ^(77)± = x{y)k ± x{y)j 



— exp 



(15) 



which is exactly the Wigner distribution (Fig. [3]). Let 
us recall that Wigner's result holds exactly for 2x2 real 
symmetric matrices, it serves as an excellent approxima- 
tion for N X N matrices though. We also know that the 
spacing distribution for a Poissonian random process in 
a plane is exactly the same as Wigner surmise. Thus 
our result proves that the complex eigenvalues of ran- 
dom cyclic matrices describe such a process. This is a 
very beautiful, non-intuitive result which brings out yet 
another characteristic of RCM. 
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following result for the JPDF for even TV: 
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FIG. 3: We observe a linear level repulsion between two 
eigenvalues which are neither real nor complex conjugate pairs 
for an ensemble of 100 x 100 matrices with 5000 realizations. 
The agreement with GOE is deceptive; in fact, this suggests 
that the eigenvalues describe a Poisson process on a plane. 

The eigenfunctions of M corresponding to the real 
eigenvalues {Ei and En_j^-^) are also simultaneously 
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eigenfunctions of "generalized parity" 77. However, the 
eigenfunctions of M corresponding to the complex con- 
jugate pair of eigenvalues are not simultaneously eigen- 
functions of Tj. Thus, when these complex eigenvalues 
occur, "generalized parity" is said to be spontaneously 
broken. Also, the eigenfunctions corresponding to the 
complex conjugate pair of eigenvalues have zero VT - 
norm. This is expected from the recent works [l|, 
on 'PT-symmetric quantum mechanics. This observa- 
tion then fully embeds our findings into the new random 
matrix theory developed recently for pseudo-Hermitian 
Hamiltonians. However, we also note that the eigenvec- 
tors ipi ('02) corresponding to complex conjugate eigen- 
values, A (A*) satisfy orthogonality defined with respect 
to rf. Since these results are found ior NxN matrices, we 
believe that this work extends the random matrix theory 
in a significant way. The findings on the spacing distri- 
butions have led us to a linear level repulsion among dis- 
tinct complex eigenvalues, whereas the spacing between 
complex-conjugate pair is Gaussian-distributed. 

Ginibre orthogonal ensemble with Gaussian dis- 



tributed real elements has been completely solved only 
recently (sEI . [sgI ] . The ensemble of asymmetric random 
cyclic matrices is a simple nontrivial instance for which 
all the interesting quantities are analytically obtained in 
an explicit manner. Such examples play an important 
role in developing a deeper insight, even when formal re- 
sults exist. 

Also, we would like to point out the role played by 
level repulsion when a system with spectral properties 
described by RMT approaches equilibrium. In its ap- 
proach to equilibrium, the central quantity of interest 
is the two-time correlation function, the long-time be- 
haviour is decided by the degree of level repulsion as the 
levels get closer. We can immediately see [13 that linear 
level repulsion is related to the exponent, 2 in t~^-tail at 
long times. 

It is a great pleasure to thank Bob Dorfman, University 
of Maryland, College Park, U.S.A. for bringing to our 
notice the role played by cyclic matrices in certain models 
in statistical mechanics, the work of Ted Berlin and Mark 
Kac [3] in particular. 
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